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We find using LSDA+U band structure calculations that the novel one-dimensional cobaltate 
Ca3Co2C>6 is not a ferromagnetic half-metal but a Mott insulator. Both the octahedral and the 
trigonal Co ions are formally trivalent, with the octahedral being in the low-spin and the trigonal in 
the high-spin state. The inclusion of the spin-orbit coupling leads to the occupation of the minority- 
spin d,2 orbital for the unusually coordinated trigonal Co, producing a giant orbital moment (1.57 
/is). It also results in an anomalously large magnetocrystalline anisotropy (of order 70 meV), 
elucidating why the magnetism is highly Ising-like. The role of the oxygen holes, carrying an 
induced magnetic moment of 0.13 fiB per oxygen, for the exchange interactions is discussed. 

PACS numbers: 71.20.-b, 71.70.-d, 75.25,+z, 75.30.Gw 



The one-dimensional (ID) cobaltate Ca3Co206 Q 
and isostructural A$MM'0§ (A=alkaline-earth metal; 
M,M'=transition metals) compounds 0] attract now 
considerable attention due to their peculiar crystal struc- 
ture and very unusual properties. A large number of re- 
search activities have been triggered by the discovery of 
the unconventional magnetic structure with strange mag- 
netization jumps [1 SI! H 0, H ©, EH El E2 , reminis- 
cent of quantum tunneling of magnetization in molecular 
magnets, and also by the observation of the large ther- 
moelectric power and magnetoresistance [l3l FLU Il5| . 

From a theoretical point of view, the understanding 
of the Ca3Co206 system is far from satisfactory. The 
crystal structure consists of triangular lattice of c-axis 
chains of alternating CoOeOctahedra and trigonal prisms 
sharing common faces jjj, |3j ■ Both the spin and valence 
state of the Co ions in these two positions, i.e. Co oc t and 
Cotrig, respectively, are a matter of controversy. There 
were several s ugg estions made Eil 
EllE2lE3lE4lE5lEflj |. the most common ones being either 
the non-magnetic low-spin (LS, S=0) Co^+ and magnetic 
high-spin (HS, S=2) Cof+ g , or the LS (S=l/2) Go 4 + t and 
HS (S=3/2) Go 2 t + g . Recent LSDA and GGA calculations 
do not provide clarity: Whangbo et al. E3 and Eyert et 
al. [lij l support the first alternative, while Vidya et al. 
E3 the second one. To add to the confusion, all those 
calculations 

E1E1E3 predict Ca 3 Co 2 06 to be a metal, 
in strong disagreement with the experiments [H ll.il Il5| . 
In fact, it has been even claimed [19| that this material 
would be the first ID oxide displaying ferromagnetic half- 
metallicity (FMHM). 

Another problem is the nature of magnetism in this 
material. In modeling the intriguing magnetic properties, 
it was always concluded or assumed that the Co chains 
behave as an Ising system with ferromagnetic (FM) in- 
trachain exchange US HE HSU [13 The nature of 
such behavior, however, is completely unclear: the ori- 
gin of Ising-like behavior was not discussed at all, and 
attempts to explain the FM exchange using an ionic su- 
perexchange model [13] meet some problems. 



In view of these controversies, we carried out a the- 
oretical study of the electronic structure and magnetic 
properties of Ca3Co206 in which we took more explic- 
itly into account the correlated motion of the electrons 
typical in transition metal oxides. In this work we used 
the LSDA+U method [2H, including also the spin-orbit 
coupling (SOC). We settle theoretically the issue on the 
valence (spin) state, and find that Ca3Co206 is a Mott or 
rather charge-transfer insulator already for moderate 
values of Hubbard U. 

The inclusion of the SOC leads to a series of surprising 
but experimentally sound results: a giant orbital moment 
for the HS Cotrig ion consistent with the experimentally 
observed large total magnetic moment, an anomalously 
large magnetocrystalline anisotropy important to explain 
the Ising-like properties, and an unusual orbital occupa- 
tion of the Cotrig ion requiring a completely new ap- 
proach for modeling the exchange interactions. 

Our calculations were performed by using the full- 
potential augmented plane waves plus local orbital 
method j^. We took the Ca3Co206 structure data de- 
termined by neutron diffraction measurements at 40 K 
P|. The muffin-tin sphere radii are chosen to be 2.7, 2.0 
and 1.6 Bohr for Ca, Co and O atoms, respectively. The 
cut-off energy of 16 Ryd is used for plane wave expansion, 
and 5x5x5 k-mesh for integrations over the Brillouin 
zone. The SOC is included by the second-variational 
method with scalar relativistic wave functions 23]. The 
easy magnetization direction is along the c-axis chains. 

We plot in Fig. 1 the LSDA results, which give a 
FMHM solution, in agreement with previous calculations 
[TR EiL Eflj | . One should note that the Fermi level is lo- 
cated in a narrow Cotrig 3d band with no more than 1 eV 
width, which is consistent with the system being ID. It is 
then also natural to expect that already modest electron 
correlation effects at the Co sites will be able to turn this 
material into a Mott insulator. 

Another important point to note from our LSDA re- 
sults is that for the Cotrig ion, the crystal fields do not 
split up the 3d orbitals into the usual e g and ti g levels. 
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FIG. 1: Density of states (DOS) of Ca3Co206 in the ferro- 
magnetic half-metallic (FMHM) state from LSDA. The Fermi 
level is set at zero energy, (a) shows (from upper to lower pan- 
els) the total DOS, the Co oct 3d, Co tr i g 3d, and O 2p partial 
DOS. (b) shows the Cotrig 3d orbitally-resolved partial DOS. 
Solid (dashed) lines depict the spin-up (down) states. 



Instead, it is found that the x 2 — y 2 is degenerate with 
the xy, and the yz with zx orbital, as can be seen from 
Fig. 1(b). In the presence of the SOC, it is then better 
to use the complex orbitals do, 0^2/^-2, an d di/d-i. Sig- 
nificant from the LSDA is also that the narrow do singlet 
and d2/d-2 doublet bands are almost degenerate. Fig. 
2 sketches the local crystal field energy diagram. This 
implies that SOC, so far not included in all the LSDA 
calculations [171 liii LUl , could have a substantial effect 
on the outcome for the predictions on the magnetic prop- 
erties. An interesting aspect from Fig. 1(b) is the finding 
that the d\jd—\ band is split off from the do and dijd-i 
bands by about 1 eV, i.e. much larger than the SOC en- 
ergy scale. Therefore, the SOC Hamiltonian can be sim- 
plified into just Q, z s z by neglecting the l+s- and Z_s+ 
mixing terms. This has far reaching consequences for the 
magnetocrystalline anisotropy as we will show below. 
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FIG. 2: Local crystal field energy diagram for: (left) Co oc 
and (right) Co tr i g without and with spin-orbit coupling. 



In Fig. 3 we show our LSDA+U+SOC results [2J. 
We find that Ca3Co20g is an insulator with ferromag- 
netic chains (FMI), i.e. not a metal, in contrast to the 
LSDA results d [H 0- We have used U=5 eV for 
the Coulomb energy and J#=0.9 eV for the Hund ex- 
change parameter. These values are taken from electron 
spectro scop y measurements [25| and Hartree-Fock calcu- 
lations 26]. To investigate to what extent the insulating 
solution is robust against the particular choice of U , we 
have calculated and plotted in Fig. 4 (top panel) the 
band gap as a function of U. One can observe that a 
value as small as 2.5 eV is already sufficient to open up 
a gap. This is consistent with the system being ID re- 
sulting in a small width of the Co 3d band. The presence 
of this gap does not practically depend on the inclusion 
of SOC and is due to the Hubbard term U. We therefore 
can conclude that Ca3Co206 is a correlated insulator j2^ 
for any realistic value of U, i.e. 4 — 6 eV for Co 3+ . 

As far as the valencies are concerned, we find from 
the LSDA+U+SOC that both the Co oct and Co tr i g sites 
have the 3+ valence. The 4+/2+ scenario (also FMI) is 
unstable and lies higher by 84 meV per formula unit (fu) . 
The Co oct ion in the 3+/3+ ground state is definitely LS: 
as can be seen from Fig. 3(a), the large Co oct t2 g -e g band 



gap does not leave any room for a 3d 6 state other than the 
one with t^g, i.e. non-magnetic. The Cotrig ion on the 
other hand, is in the HS state (S=2). Very interesting is 
the orbital occupation for this Co^ ion: Fig. 3(b) shows 
that five electrons occupy all the spin- up do, di/d-i, and 
d2/d_2 orbitals, and that the sixth electron resides in the 
spin-down di orbital. This is the consequence of the SOC 
in this very narrow band system with the tiny crystal field 
splitting between the do and d^jd-i levels. When the d_2 
orbital is forced to be occupied instead, the increasing 
total-energy allows us to estimate the SOC strength to be 
about 70 meV. As a result, the orbital contribution to the 
magnetic moment can be much larger than the 1 [Ib one 
usually finds for Co or Fe ions in Oh symmetry. An ionic 
d2 state would give an orbital moment of 2 /x^, but here 
we find 1.57 hb due to covalency and small distortions 
of the trigonal prism, see also Table I. Experimentally, 
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FIG. 3: Density of states (DOS) of Ca 3 Co 2 06 in the ferro- 
magnetic insulating (FMI) state from LSDA+U+SOC. 



Hu et al. have recently found using soft-x-ray magnetic 
circular dichroism that Ca3Co20g indeed has an orbital 
moment which is significantly larger than 1 /is [27^ . 

Adding the orbital moment to the 4 [Ib from the S=2 
spin contribution, we end up with the total magnetic 
moment of about 5.66 /x^/fu. Experimentally, Maignan 
et al. 9] found from magnetization measurements the 
saturated total moment close to 5.0 (j,b- 

We may speculate that perhaps a small part of the 
sample («10%) is somewhat misaligned, which due to 
the anomalously large magnetocrystalline anisotropy as 
we will show below, will then not get magnetized. In any 
case, since it is easier to find experimental reasons why 
the sample is not fully magnetized rather than 'over'- 
magnetized, it is probably better that the theory slightly 
overestimates the moment. It would be especially un- 
satisfactory, if one, for example, would ignore the SOC 
or assume that the sixth electron occupies the do spin- 
down orbital |28| . giving in both cases a total moment 
of only 4.0 /j,b- In fact, as listed in Table I, we find in 
our LSDA+U+SOC calculations that the configuration 
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FIG. 4: U dependence of the insulating gap (top panel), as 
well as of the spin and orbital contributions to the Cotrig 
magnetic moment (bottom panel). 



with the sixth electron forced to occupy the do spin-down 
orbital is unstable by 71 meV/fu as compared to the 
ground state. Important is that this analysis is also ro- 
bust against the particular choice of U made here. Fig. 4 
(bottom panel) shows that the spin and orbital moments 
stay constant within 0.2 \xb as long as the system is an 
insulator, i.e. when U is varied between 2.5 and 6 eV. 

The occupation of the spin-down c?2 orbital has as a 
consequence that the orbital and, due to the SOC, also 
the spin contributions to the magnetic moment are ori- 
ented along the z-direction, i.e. the c-axis. To flip both 
the spin and orbital moments into the plane perpendic- 
ular to the c-axis is impossible, since this would require 
that the spin-down d\ and d-\ become partially occu- 
pied, and we have seen from the LSDA calculations in 
Fig. 1(b) that this would cost roughly the crystal field 
energy (~1 eV). To flip only the spin moment into the 
plane but keep the orbital moment along the c-axis is also 
quite difficult, since this would cost the full spin-orbit 
splitting of about 70 meV. Also a third scenario in which 
one tries to alter the moment directions by occupying 
non-spin-orbit-active spin-down do orbital is equally un- 
likely, since our LSDA+U+SOC calculations show that 
this would cost 71 meV, as already mentioned above and 
in Table I. All this means that the magnetocrystalline 
anisotropy is exceptionally large, and that all the rel- 
evant magnetic degrees of freedom are highly fixed in 
the z-direction, justifying why the magnetic behavior of 
Ca3Co206 can be well described by Ising models. 

In trying to explain the intra-chain FM exchange in- 
teraction we also have investigated the antifcrromagnetic 
insulating (AFI) scenario. Here we took the chain as an 
alternation of spin-up and spin-down HS-CotH g ions with 
LS-Co oct ions in between. Our LSDA+U+SOC calcula- 
tions find that the AFI solution, having the same spin 
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TABLE I: Calculated electronic state of CasCo206; total energy difference (meV) per formula unit; moments (/j-b) at each 
Cooct, Cotrig, and O ion, as well as in the interstitial region and total magnetic moment per formula unit. Both Co sites are 
trivalent in the solutions listed. The high-spin Co|^ fl has either the minority-spin d,2 or do occupied. 
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and orbital moments at the HS Cotrig ions as the FMI 
ground-state, lies above the latter by about 12 meV/fu, 
see Table I 29]. As a result, the intra-chain FM coupling 
parameter of each HS Cotrig pair can be estimated to 
be of the order of 1.5 meV (17 K), assuming a simple 
Heisenberg model with S=2. This intra-chain exchange 
parameter is in reasonable agreement with the experi- 
mentallyobserved intra-chain Curie temperature of 24±2 

kE!1. 

This result seems at first sight to support the ionic 
superexchange model proposed by Fresard et al. [2fj| . 
That model used an ansatz for the Co ions which was 
confirmed by our LSDA+U+SOC, namely LS Co 3 + t and 
HS Co 3 + ff with the do orbital half occupied, and found 
the exchange parameter of about 2 meV, very close to 
ours. However, that model considered basically a set of 
virtual excitations which are purely ionic in character and 
involve essentially only the Cotrig and Co oc t 3z 2 — r 2 or- 
bitals, in which case the solution should be AF, according 
to the first Goodenough-Kanamori- Anderson rule [3(J • 

We believe that the intra-chain FM interaction is con- 
nected with the large contribution of holes on the oxy- 
gens. On the basis of general arguments and also our 
calculations, one finds that the actual charge distribu- 
tion corresponds rather not to Co 3+ , but to Co 2+ L 
(L=oxygen hole), although the quantum numbers are 
those of Co 3+ . In this case two of such oxygen holes 
from two neighbouring HS-Co 3 ^ g form a triplet state at 
the LS-Co^ in between, much the same as two oxygen 
holes form a triplet state around a metal vacancy in ZnO 
and CaO (U- This allows to gain the full coherent hy- 
bridization energy with the empty e g states of the Co 3 |t t 
ion. An indication for this is the presence of a covalency- 
induced spin moment in the oxygen of about 0.13 \xb 
per ion as well as in the LS-Co 3 ^ of about 0.07 fj,s, see 
Table I. This picture of the origin of the FM intra-chain 
exchange should be confirmed by further studies. 

To conclude, we find using LSDA+U+SOC band struc- 
ture calculations that Ca3Co2 0g is not a ferromagnetic 
half-metal but a Mott insulator with both the octahe- 
dral and the trigonal Co ions being formally trivalent, 
which settles the valence-state issue. Spin-orbit coupling 
and unusual coordination of the trigonal Co ion lead to 
the occupation of the d,2 spin-down orbital, generating a 



giant orbital moment (1.57 hb) and an extremely large 
magnetocrystalline anisotropy (70 meV), which explains 
the Ising character of the magnetism of this material. 
The ferromagnetic intra-chain interaction is presumably 
connected with the effect of oxygen holes in this small 
charge-transfer gap system. This also leads to an appre- 
ciable magnetic moment of 0.13 \xb on each oxygen. 
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